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Abstract. Smooth real cubic surfaces are birationally trivial (over R) if and only if their 
real locus is connected or, equivalently, if and only if they have two skew real lines or two 
skew complex conjugate lines. In such a case a parametrization over the reals can be given by 
cubic polynomials. In this short note we provide a simple geometric method to obtain such 
parametrization based in an algorithm by I. Polo-Bianco and J. Top [14] . 



Introduction 

A cubic surface S is the vanishing set of a homogeneous polynomial / of degree 3 in P 3 , i.e., 

S = {{x : y : z : t) G P 3 | f{x : y : z : w) = 0}. 

By a classical result of Clebsch [7] we know that a smooth cubic surface over the complex 
numbers admits a parametrization by cubic polynomials defined over C (see Theorem [2] below). 
If S is a cubic surface defined over the real numbers, then S admits a parametrization defined 
over the real numbers if and only if the real locus of S is connected or, equivalently, if and only 
if S has two disjoint real lines or two disjoint complex conjugate lines (see Theorem [3] below) . 

Algorithms for parametrizing smooth cubic surfaces have shown to be of great interest in 
Computer Aided Geometric Design since these surfaces are much more flexible than quadrics 
and the resolution of many problems in science and technology depends on using the so called bi- 
cubic patches [T] (bounded and nonsingular cubic surfaces in M 3 ) to approximate certain objects 
(such as molecules, organs of the human body, etc.). Some applications of these particular cubic 
surfaces (A-cubic patches) to molecular modeling can be found, for example, in [3l Hj. Some 
of the known algorithms for parametrizing cubic surfaces can be found in the works of T. W. 
Sederberg and J. P. Snively [15], C. Bajaj, R. Holt and A. Netravali [2j, T. G. Berry and R. R. 
Patterson [6] or I. Polo-Bianco and J. Top [H]. 

In this note we parametrize smooth cubic surfaces by giving a purely geometric and very 
simple construction based on the algorithm presented in [14J. The main advantage of this new 
approach is that it is easy to understand and easy to program, as we will show at the end 
of this note. This a is a small, but we believe important contribution to the theory of real 
cubic surfaces, that might interest the mathematical community in general and in particular 
the CAGD researchers. The construction presented here is, from the computational point of 
view, optimal. It requires the computation of the lines on a smooth cubic surface, a problem 
recently solved by R. Pannenkoek in his master thesis (compare [14, §3.1] and [I2])- However it 
does not require a different construction for the cases where the cubic surface has two disjoint 
real lines or two disjoint complex conjugate lines, as the previous algorithms did (see [2], [6], 
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[14J). We have implemented the algorithm in the mathematical software MAPLE and applied 
it to the Fermat cubic surface and to a generic smooth cubic surface (see Section HJ). 

1. Cubic surfaces 

In 1949 A. Cayley and G. Salmon discovered that smooth cubic surfaces over C have precisely 
27 lines over C. The 27 lines in a smooth cubic surface have a very special symmetry and 
configuration and whole books have been written to describe them (e.g. [8], [IS]). L. Schlafli 
introduced the concept of a double six in order to describe the intersection behaviour of the 27 
lines in a smooth cubic surface. A double six is a set of 12 of the 27 lines on a cubic surface S, 
represented with Schlafli's notation as: 

/ a\ ci2 0,3 04 a 5 a§ \ 
^ 61 b 2 63 h h b 6 J ' 

with the following intersection behaviour: a* (resp. fej) does not intersect a,j (resp. bj) for % ^ j 
and cii intersects bj if and only if % ^ j. The plane spanned by two intersecting lines and bj 
for i ^ j intersects the cubic surface on a third line Cy. Repeating this process the remaining 
15 lines on S are obtained together with their intersecting behaviour. 

The following property regarding the lines on S will be needed in what follows, and can be 
easily verified by using the configuration of the 27 lines on S. 

Lemma 1. Let S be a nonsingular cubic surface and £\ and £2 are two skew lines on S. Then 
there are precisely 5 lines that intersect both l\ and £2- 

Proof. See [HJ Lemma 1.2] □ 
Another important result concerning smooth cubic surfaces was obtained by A. Clebsch in 

Theorem 2 (Clebsch). Let S C P 3 (C) be a nonsingular cubic surface. Then S is obtained by 
mapping P 2 to P 3 by the space of cubic forms passing through 6 points in general position (i.e. 
no three on a line and no six on a conic). 

Proof. See [3 Proposition IV. 12]. □ 
A base change is just a projective transformation of the cubic surface. 

2. Real cubic surfaces 

If S is a real cubic surface then it is not true in general that the morphism in Theorem [2] can 
be defined over the reals. For example, the cubic surface defined by the following equation is 
not parametrizable since its real locus is not connected (see [T3]). 

(1) X 3 + Y 3 + Z 3 + X 2 Y + X 2 Z + XY 2 + Y 2 Z - QY 2 W + 11ZW 2 - 6W 3 = 0. 

One has the following characterization of birationally trivial real cubic surfaces. 

Theorem 3. Let S be a smooth cubic surface defined over R. Then the following conditions 
are equivalent 

(1) The real locus S{R) is connected. 

(2) S has a set of two skew lines defined over R (i.e. either both real or complex conjugate). 

(3) S is birationally trivial over R (i.e. S is parametrizable over R). 
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(4) S is the blow up of P 2 at six points denned over R (i.e. S admits a parametrization by 
cubic polynomials). 

Proof. Compare pi], [ID], [9] and p3]. □ 



3. The parametrization 

The geometric construction we propose works as follows. Let 5* be a smooth cubic surface 
defined over C and 4 and 4 two skew lines on S (i.e. 4 fl 4 = 0)- Consider m a line on S that 
intersects both 4 and 4 and H a plane containing m that does not contain 4 or l 2 . Then for 
any point x G H there exists a unique line i x passing through x and intersecting both 4 and 
£2- The line £ x intersects the smooth cubic surface at three points: namely £ x H4, 4 H4, and 
a third point that we will denote by q x . Now, the parametrization is given by the map 

defined by $(x) = q x . 

A real smooth cubic surface S admits a parametrization defined over R if and only if there 
exists a pair of skew lines on S (call them £\ and £2) both real or complex conjugated (see The- 
orem [3]). Let us show that, for such S, the parametrization $ defined above can be constructed 
over the real numbers. Indeed, there are precisely five lines on S intersecting both 4 and £2 
(see Proposition [T]) where at least one of them, call it m, is real. This implies that the plane H 
can be chosen to be real, and therefore, the map $ will be defined over R. 

We continue with the proof that the construction proposed above gives a parametrization of 
5 by cubic polynomials. The proof is purely algebraic and it is not needed for the application 
of the geometric construction to concrete examples. 

Theorem 4. Following the above notation, the map $ : P 2 — > S C P 3 is a birational morphism 
given by cubic polynomials. 

Proof. The morphism $ factor as follows: 

where $i(a;) = {£ x D 4, 4 H £ 2 ) G 4 x 4 and $ 2 (4 H 4,4 H 4) = <?x- The morphism $1 blows 
up the points £±H H and 4 H iZ and blows down the line m. The morphism (I?^ 1 blows down 
the five lines that intersect both £1 and 4, in particular the line m. Therefore $ = $ 2 
blows up i7 at six points. Hence $ is given by cubic polynomials. □ 

Remark 1. In the last section of this text we give a computational proof of this fact using the 
program MAPLE. 

Remark 2. (a) If both lines i\ and 4 are defined over R and one considers the map $2 : 
4x4 = P 1 x P 1 — > S in the proof of the previous theorem, then $2 is a parametrization 
of S by biquadratic polynomials. In fact $ 2 is defined by the space of bihomogeneus 
forms of bidegree (2, 2) passing through 5 points qi , . . . , q§ in general possition in P 1 x P 1 . 
This is the algorithm proposed in [T5] and [2] with a slight modification for the case 
when S contains two complex conjugate disjoint lines. The five points qi, . . . ,q 5 are the 
base points of the parametrization. 



4 



JON GONZALEZ-SANCHEZ AND IRENE POLO BLANCO 



(b) If the plane H in the construction above does not contain a line in S, then the map $ 
is given by polynomials of degree 5. In fact, in such a case $ is defined by the space 
of forms of degree 5 vanishing with multipicity two at five points q\, . . . ,q§ and with 
multiplicity one at two points p\ and p 2 such that q±, . . . , q^, p\ and P2 are in general 
possition. The seven points qi, . . . , q5,Pi,P2 are the base points of the parametrization. 



4. Examples 

We will now apply the construction given in the previous section to some concrete examples. 
We start with the Fermat cubic. 

4.1. Fermat cubic. The Fermat cubic surface S in M 3 is defined by 1 + x\ + x\ + x\ = 0. 
Following the notation of the previous section we take £\ = (— uj 2 , —tot, t) and £ 2 = (-co, —u 2 t, t) 
where a; is a primitive third root of unity [HI §5.2]. The line m is (— t, — l,t) and we choose 
the plane H to be x 2 = — 1. Now the algorithm that parametrizes S can be easily implemented 
in a computer algebra program (e.g. in MAPLE). 

# Select the Fermat cubic: 

f:=(xl,x2,x3)^ l+xr3+x2~3+x3~3: 
w:=RootOf(x"2+x+l,x): 

# Choose the lines £\ and £ 2 on S: 
linel :=r — > [-w"2, -w*r, r]: 

Iine2 := s — > [-w, -w"2*s, s]: 

# The line £ x : 

linex := t [yl,-l,y2] + [t, t*b,t*c]: 

Sol := solve({linel(r)[l] = linex(tl)[l], linel(r)[2] = linex(tl)[2], linel(r)[3] = linex(tl)[3], 
Iine2(s)[l] = Iinex(t2)[l], Iine2(s)[2] = Iinex(t2)[2], Iine2(s)[3] = Iinex(t2)[3], f(linex(t3)[l], 
Iinex(t3)[2], Iinex(t3)[3]) = 0, tl <> t3, t2 <> t3}, [r, s, b, c, tl, t2, t3]): 

# The pametrization of S is 
parametrization := factor(subs(Sol[l], Iinex(t3))); 

This program gives the following parametrization of S: 

yf - iy\ - ybi + 3?/i + Zym + y\y\ + y\ 



xi(yi,y2) 
^2(^/1,1/2) 
^3(2/1, 2/2) 



yf - 2 2/i - 2/i2/2 + 3yi + 2y x y 2 + y\y x - 3 - 3y 2 - 2y\ ' 
gj/i + y\y<i - yim - 2ym -3y 1 + 3 + 3y 2 + yj + 2y\ 

y\ - 2 yf - y\yi + 3 ?/i + 2 2/i?/2 + y\y\ - 3 - 3y 2 - iy\ ' 
vlvi - 2/1 - 2ym - yhx + yl + 2yf + 3y 2 

y\ - 2 2/i 2 - 2/i 2 2/2 + 3yi + 2yiy2 + 2/I2/1 - 3 - 3y 2 - 1y\ ' 



4.2. A general construction with MAPLE. Following Section[3]we consider a triple (£i,t 2 , m) 
where £\, £ 2 and m are three lines in P 3 such that £\ and £ 2 are skew and m intersects both £\ 
and £ 2 . We say that a cubic surface admits a triple (£1,^2,^) if the three lines £1, £ 2 , m are 
contained in 5*. If K is a subfield of C we say that the triple (£1,^2,^) is defined over K if 
£\ U£ 2 Um is defined over K, i.e. m is defined over K and l\ and £ 2 are conjugate on a quadratic 
extension L of K. 

For our purposes it is enough to consider the case m = (t, 0, 0), £\ = (a + bit, b 2 t, b 3 t) and 
£2 — (c + dit,d 2 t,dst) (which is the general case up to a projective transformation). First we 



BIRATIONALLY TRIVIAL REAL SMOOTH CUBIC SURFACES 



5 



compute the space V of cubic surfaces admitting the triple (£i,£ 2 ,m). If the triple (£i,£ 2 ,m) 
is defined over K so will the space V. We can easily compute the space V . 

with(PolynomialTools): # Consider a Generic cubic: 

f:=(xl,x2,x3)^A*xl^3+B*x2^3+C*x3^3+D*xl^2*x2+E*xl /s 2*x3+F*x2 /v 2*xl 
+G*x2~2*x3+H*x3~2* xl+ J*x3~2*x2+K*xl*x2*x3+L*xr2+M*x2~2+N*x3~2+ 
0*xl*x2+P*xl*x3+Q*x2*x3+R*xl+S*x2+T*x3+U; 

# Choose two disjoint lines l\ and £ 2 . 
Iinel:=t^[a+bl*t,b2*t,b3*t]; 
Iine2:=t^[c+dl*t,d2*t,d3*tj; 

# Choose a line m that intersect both l x and £ 2 : 
linem:=t-i>[t,0,0]; 

# Compute the Coefficient Vector of £ x , £ 2 and m: 

El:=expand(subs(xl=linel(t)[l],x2=linel(t)[2],x3=linel(t)[3],f(xl,x2,x3))); 

E2:=expand(subs(xl=line2(t)[l],x2=line2(t)[2],x3=line2(t)[3],f(xl,x2 

Em:=expand(subs(xl=linem(t)[l],x2=linem(t)[2],x3=linem(t)[3],f(xl,x2,x3))); 

Vl:=CoefficientVector(collect(El,t),t); 

V2:=CoefficientVector(collect(E2,t),t); 

Vm :=Coeff icientVector(collect( Em ,t) ,t) ; 

# Compute the variety of the cubics passing through the lines £\, £ 2 and m. 
Sols:=solve(Vl[l]=0,Vl[2]=0,Vl[3]=0,Vl[4]=0,V2[l]=0,V2[2]=0,V2[3]=0,V2[4]=0, 
Vm[l]=0, Vm[2]=0,Vm[3]=0,Vm[4]=0,A,B,C,D,E,F,G,H,J,K,L,M,N,O,P,Q,R,S,T,U): 
Cubic: =subs(Sols,f(xl,x2,x3)); 

This program gives a description of the space of cubics passing through the lines £±, £ 2 and 
m. One can proceed by computing a parametrization of a general cubic in V with the method 
proposed in Section [3] as follows. 

# Consider the plane H=(yl,yl+y2,l-y2) 
linex := t [yl,y2,0] + [t, n*t, m*t]; 

# Solve 

Sol := solve(linel(r)[l] = linex(tl)[l], linel(r)[2] = linex(tl)[2], linel(r)[3] = linex(tl)[3], 
Iine2(s)[l] = Iinex(t2)[l], Iine2(s)[2] = Iinex(t2)[2], Iine2(s)[3] = Iinex(t2)[3], subs(xl=linex(t3)[l], 
x2=linex(t3)[2], x3=linex(t3)[3], Cubic)=0, tlot3,t2ot3, [r, s, n, m, tl, t2, t3]): 

# Now the parametrization is 
parametrization := factor(subs(Sol[l], Iinex(t3))): 

This gives a parametrization by cubic polynomials (in y\ and y 2 ) as it is claimed in Theorem 

H 



^i (2/1, 2/2) 
2:2(2/1,2/2) 
^3(2/1,2/2) 



-2Ta 2 bldld 2 yl - TVgcMg^yg - 3Ta 2 b 2 cd 3 b 3 d 2 yi + ... 
2Ca& 2 cd|&i6 3 y| + 2Ka 2 b 2 cb\d\y\y 2 — 2Ka 2 b 2 cd\h\d iy \ + . . . ' 
2d 2 b\cTb x d\y x yl - 2d 2 b\acd\Pb x y x y\ - 2d 2 b\aTb 1 d\y x yl + ... 
2Cab 2 cd\%\b 3 y\ + 2Ka 2 b 2 cb\d\y\y 2 — 2Ka 2 b 2 cdl 2 d)\d\y\ + . . . 
b 3 d 3 ab\(?Ddxd 2 y\ — 3b 3 d 3 a 2 b 2 c 2 Ddid 2 y 2 + 2b 3 d 3 ab\e'Dd\yxy 2 + . . . 
2Cab 2 cd 3 bib 3 y 2 + 2Ka 2 b 2 cb\d\y\y 2 — 2Ka 2 b 2 cdl 2 d)\d\y\ + . . . 
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